Force-Extension Relations for Polymers with Sliding Links 
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Topological entanglements in polymers are mimicked by sliding rings (slip-links) which enforce 
pair contacts between monomers. We study the force-extension curve for linear polymers in which 
slip-links create additional loops of variable size. For a single loop in a phantom chain, we obtain 
exact expressions for the average end-to-end separation: The linear response to a small force is 
related to the properties of the unstressed chain, while for a large force the polymer backbone can 
be treated as a sequence of Pincus-de Gennes blobs, the constraint effecting only a single blob. 
Generalizing this picture, scaling arguments are used to include self-avoiding effects. 
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Entanglements play an important role in the behav- 
ior of macromolecules. For instance, mechanical links 
(e.g., in catenanes) and knots naturally appear in long 
polymers [Q. In biological systems, specific proteins act 
upon topological states: the degree of entanglement of 
chromosomes during cell division or knotted states 
in bacterial DNA which may arise during random ring 
closure [9, can be modified by topological enzymes jjj. 
Synthetic RNA trefoil knots have been used to prove 
the existence of a similar, previously unknown, topology- 
changing enzyme J5j. Tight molecular knots have even 
been found deep inside the native state of proteins 
Experimental advances now make it possible to manip- 
ulate single molecules by optical tweezers. Thus, tight 
knots could be tied into single actin filaments or DNA 
strands 0. Mechanical properties, and forces in the pN 
range relevant to biopolymeric processes, can be mea- 
sured by atomic force microscopy, or more direct mi- 
cromechanical methods [||. It is therefore possible to 
record the force-extension (FE) curve of single polymers 
with a fixed topology, from which valuable information 
about the properties of a molecule can be obtained and 
compared to theoretical predictions. 

While there has been extensive progress in the statis- 
tical mechanics of polymers in the last decades @, [To| , 
the analysis of topological constraints is hampered by 
the difficulty of treating the resulting division of phase 
space into accessible and inaccessible regions. Since the 
mathematical methods of knot detection using topolog- 
ical invariants JHJ cannot be conveniently incorporated 
into a statistical-mechanical formulation, one may try to 
use geometrical constrictions to mimic knots: Consider 
a polymer threaded through a small ring as depicted in 
Fig. pi and not allowed to withdraw from it, although the 
ring may freely slide along the polymer and the loop size 
can change. Constrictions of this type (called slip-links 



FIG. 1: Polymer threaded through a slip-link (ring) forms 
a loop and two dangling ends. It is not allowed to withdraw 
from the slip-link. 



(SLsV) were introduced by Ball, Doi, Edwards and oth- 
ers |Q to investigate the elasticity of rubber, where they 
were used to represent entanglements between different 
polymers. Recently, a detailed study of the size distri- 
bution of loops in SL structures was performed Jl3| . In 
this work we consider FE curves of such polymers, with 
or without self-avoidance. We show that the knowledge 
of the statistics in the absence of the force, combined 
with the Pincus-de Gennes ]l3] blob picture, suffices to 
understand many features of the FE curves. 

We first discuss the statistical effects of the slip-link for 
non-self- avoiding chains (phantom polymers). Consider 
an A-stcp chain with mean-squared step size a 2 in d- 
dimensional space. In the N ^S> 1 limit, the probability 
density function for the end-to-end distance (EED) r is 
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This expression also describes an A-step random walk 
on a (hyper)cubic lattice with lattice constant a. The 
number of closed TV-step loops on such a lattice is 
{2d) N [d/(2nN)] d / 2 . Now consider the case where a 
SL forces a closed loop of size n. Since the loop 
can be located on any of N — n positions on the TV- 
step chain, the number of possible configurations of 
the combined system of the polymer and the SL is 
(A - n)(2d) N - n (2d) n [d/{2Tm)} d / 2 . Thus, for n > 1 and 
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(N — n) 3> 1, the probability for a given n is 
3Ar (n) = Af{N - n)n- d/2 , 



(2) 



where Af is a normalization factor. For d > 2, Af depends 
on the short-distance cut-off, i.e., on the microscopic de- 
tails of the walk |l5| . The presence of the SL also modifies 
the probability density of the EED, p(r), by reducing the 
length of the backbone to (N — n) . 

When the end-points of a polymer are stretched by a 
force f , its properties can be derived from the partition 
function 



J d d vp(v) e f r / T , 



(3) 



where / = |f|, and the temperature T is in energy units, 
i.e., kg = 1- In the presence of the force, the vector r 
is on average parallel to f, and from Eq. (||), its mean 
projection along f is given by 



(r)f = T 



d\nZ 



(4) 



In particular, for p{v) = pn(j) in Eq. ([!]), the parti- 
tion function is Z = exp[iVa 2 / 2 / (2dT 2 )], while the mean 
EED is a linear function of /, (r) f = fNa 2 /(dT), for 
any arbitrary value of the force. Since the mean-squared 
EED R 2 of an unstrained (zero-force) phantom chain is 
Na 2 , the FE relation can be re- written as 



' )f = dT f - 



(5) 



More complicated forms of p(r) do not lead to a simple 
linear relation, and in many cases the relation between 
(r) f and an arbitrary / cannot be calculated exactly. 
However, for sufficiently small /, linear response theory 
provides a simple universal answer: By expanding the ex- 
ponent in Eq. (j^) in powers of /, and by omitting powers 
higher than 2, we see that Eq. ([|) is valid for arbitrary 
spherically symmetric p(r) , provided that R 2 is the mean 
squared EED calculated at zero force. The force can be 
considered small when (r) f <C R, i.e., for / -CT/i?. 

The probability density of the EED of a phantom chain 
with a simple SL (Fig. [l]) is given by 

p(r) = PN-n{r)gN{n) (6) 

(see Eqs. © and (§)). Here, N — n is the number of 
monomers in the force- carrying backbone of the polymer. 
Thus, R 2 can be found by integrating r 2 with the above 
statistical weight, over all possible r and n, leading to 

R 2 =a 2 (N- 
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FIG. 2: Force-extension curves for a 100-monomer polymer 
with a slip-link in d — 3 with its loop threaded through addi- 
tional (from top to bottom) m — 0, 1, 2, 3 and 4, rings. 



Here, Cd are short length-scale cutoff-dependent con- 
stants. In this expression, (n)o is simply the mean num- 
ber of monomers inside the constricted loop, and R 2 is 
obtained by replacing N in the expression for R 2 of a 
simple phantom chain by N — (n)o- We note that for 
d > 2, (n)o has a sub-linear dependence on N, and as d 
increases the correction created by the SL depends more 
weakly on N JL6| . For d > 4 random walks do not form 
long loops and (n)o becomes independent of N. 

The expression for R 2 of a phantom chain with a SL 
can now be substituted in Eq. (^), to obtain the FE re- 
lation (r) f — f(a 2 /dT)(N — (n)o) for small /. For large 
/, this expression is no longer valid. However, by direct 
inspection of the required average we find that Eqs. (|), 
(§), and (§) lead to 
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where (n) / is the mean loop size in the presence of the 
force, equal to 



dnng N (n) exp [- f 2 a 2 n / {2dT 2 )} 



dng N (n) exp [-Pa*n/(2dT*)] 



(9) 



The lower limit of no in the above integrals is the minimal 
loop size allowed by the specific model. 

The FE relation is no longer linear, even for a phan- 
tom polymer with one SL, although the deviation dis- 
appears rapidly with increasing /. For the simple SL 
with the weight in Eq. (||), the nonlinearity is barely de- 
tectable even for the relatively small value of N = 100, 
as indicated by the top line in Fig. ||. In more compli- 
cated topologies of many SLs, we frequently encounter 
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FIG. 3: A slip-link whose loop is threaded through an addi- 
tional m — 3 rings which are restricted to stay on the loop. 



the behavior <?Ar(n) ~ n~ a |L3|, where n is the total 
number of monomers that do not belong to the direct 
path between the ends of the polymer, while N appears 
in the prefactor or other non-singular parts of the prob- 
ability density. In such cases, the non-linearities become 
more pronounced as a decreases. Let us specifically con- 
sider the "toy-example" of a SL in which m additional 
rings slide around the loop, as depicted in Fig. [| The 
number of ways of placing these sliding rings leads to 
gj\i(N) — Af(N - n)n m ~ d ' 2 0. Depending on the val- 
ues of m and d, three different behaviors can be distin- 
guished: 

(i) For m > d/2—1, the integrals in Eq. (||) are dominated 
by large n. We can thus set the lower limit of the integrals 
to 0, and introduce the new variable x = n/N , to get 



(r) f _ faVN 



Jo dx '< 
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This FE curve now satisfies a scaling form (r) f/R = 
<f> ( f R/[\/2dT]\ , where the scaling function has the lim- 



its 



2/(1 + m- d/2), z^O, 
z — c/z, z — > oo. 



(11) 



The initial slope is reduced for larger m, as depicted in 
Fig. ^, while the asymptotic form at large / is reached 
with a correction that falls off as 1//. The physical origin 
of the nonlinearity is the tightening of the initially large 
loop in the intermediate regime. 

(ii) For d/2 — 2 > m > d/2 — 1, the mean loop size 
grows as 7V 2+m ~ d / 2 , playing the role of an additional 
(sub-leading) length scale. Consequently, the FE curve 
no longer has a simple scaling form, and behaves as 
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This case may be most similar to that of knots in three 
dimensions, and indeed such corrections to scaling were 
used in Ref. flUI to extract the size of the knot. 




FIG. 4: In the Pincus-de Gennes scenario, the stretched poly- 
mer is viewed as a linear sequence of "blobs" (circumscribed 
by circles). Within each blob the polymer is unstressed; the 
size and number of blobs depends on the stretching force. 



(Hi) For m < d/2 — 2, both integrals in Eq. (^) are domi- 
nated by the short distance cutoff, resulting in (n) f ss no, 
independent of /. The FE curve is thus linear in this 
regime, with finite size corrections that disappear as 1/N. 

While the above results are easily obtained for a phan- 
tom polymer with a SL, it is convenient to restate them 
in a form that is more generally valid, and, in particular, 
applicable to interacting polymers. This will be done us- 
ing the Pincus-de Gennes picture [fl4| , according to which 
a stretched polymer (without a SL) at short scales does 
not feel the influence of the external force, and correla- 
tions remain as in the unforced polymer, while at longer 
distances it is essentially a linear object aligned to the 
force. The polymer can then be visualized as a linear 
chain of blobs, as depicted in Fig. ||. The number of 
monomers Nf, inside a blob is determined by the condi- 
tion /i?^ « T, where Rb is the EED of Nb monomers. In 

1 /2 

the case of a phantom polymer (Rt = aN b ' ), this leads 
to N\j w (T/fa) 2 , while in the more general case with 
R b = aN", we get N b « (T/fa) 1 ^. Consequently, for 
large forces the EED of a whole polymer is the size of a 
single blob times the number of blobs, i.e., 



(r) f = {N/N b )Rb = aNN"- 1 = aN(fa/T)~- 



(13) 



For phantom polymers the FE curve remains linear even 
for large /, while for a self- avoiding polymer in d = 3 
(with v « 0.58) the relation is highly non-linear. 

We now note that Nb is the scale over which the ex- 
ponential factor in Eq. (^) decays, and consequently the 
mean size of the SL loop for / 3> T/ (ayN) can be esti- 
mated as 



J no b dnng N (n) 



(»>/ 



J^dng N (n) 



(14) 



However, this expression is exactly the size of the link in 
a polymer consisting of Nb monomers in the absence of 
an external force, i.e. 



(n)f,N ~ (n)o,N b 



(15) 



The first subscript in this equation denotes the size of 
the force, while the second index indicates the total num- 
ber of monomers. We can, therefore, view the SL loop 
as being confined to a single blob. Since within a blob 
the external force is not felt, its size is determined by 
regarding the entire polymer length as Nb, as depicted 
qualitatively in Fig. [| 
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FIG. 5: Qualitative representation of a stretched polymer 
with a SL. The loop created by the SL is contained within a 
single blob, and its size is determined only by the number of 
the monomers in the blob. 

While the leading term in the expression for the EED 
of a strongly stretched polymer with a SL will still have a 
form given by Eq. ( [l3l) , subleading corrections depend on 
the influence that the presence of SL has on the EED of 
the unstretched polymer. For a phantom chain there is 
a clear separation between the segment that creates the 
link and the remainder of the chain, and, consequently, 
in the absence of an external force, the reduction in R 2 
can be simply related to the reduction of N by Eq. (0). 
This does not have to be the case in the presence of inter- 
actions, and each case must be considered separately. In 
the presence of self-avoiding interactions (for d < 4), it 
can be shown jl3| , that for n <C N the statistical weight 
of an n-monomer loop is given by gN(n) ~ n ~ c , where 



c > 2. Consequently, (n)o is independent of N, and will 
cause no detectable modification in the FE curve of a 
self-avoiding polymer. This picture can be easily general- 
ized to a sequence of (non-interpenetrating) SLs. At zero 
force, the loops will compete for the available length, each 
acquiring a fraction of the overall length, as described in 
Rcf. |ll|. At strong force, in the blob regime, each slip 
link is most likely confined to its own blob. 

Viewing the stretched polymer as a sequence of blobs, 
with only individual blobs affected by the presence of the 
constraints such as SLs, creates a convenient framework 
for evaluation of FE relations. This picture may possi- 
bly be extended to knotted polymers: If the size of each 
prime knot factor depends on the number of monomers 
N as a power law N f , then the application of a strong 
stretching force will confine the knot to a blob and will 
reduce its size to iVjf. Such a scenario (at small forces) 
was explored in a recent Monte Carlo study fl§|| . 
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